AD-A148  128  AN  EXACT  SOLUTION  TO  THE  THO  POINT  OC  1/1 

(OPERATING-CHARACTER ISTIC)  CURVE  DE  .  (U)  HAROLD  LEONARD 
STUART  SCHOOL  OF  MANAGEMENT  AND  FINANCE  CHICAG. 

UNCLASSIFIED  J  E  HALL  ET  AL.  NOV  83  N00014-82-K-0560  F/G  12/1  NL 


AN  EXACT  SOLUTION  TO  THE  TWO 
POINT  OC  CURVE  DETERMINATION 

James  E.  Hall* 

Thomas  w.  Knowles  ** 

M.  Zia  Hassan  ** 


November,  1983 


AN  EXACT  SOLUTION  TO  THE  TWO 
POINT  OC  CURVE  DETERMINATION 


Janes  E.  Hall* 
Thomas  W.  Knowles  ** 
M.  Zia  Hassan  ** 

November,  1983 


Stuart  School  Working  Paper  No.  8301 


*University  of  Illinois  at  Chicago 
♦•Illinois  Institute  of  Technology 


This  document  has  been  approved 
for  public  release  and  sale;  u* 
distribute  j»  unlimited. 


ABSTRACT 


\ 

The  determination  of  an  operating-characteristic  (OC)  curve  which 
passes  through  two  arbitrarily  selected  points  Is  a  well  known  problem 
in  statistical  quality  control.  Ordinarily  this  problem  cannot  be 
solved  exactly  as  the  solution  sampling  plan  must  be  defined  by  a  non¬ 
negative  Integer  couple  (n,c). 

In  this  paper  a  randomization  process  Is  defined  for  the  acceptance 
and  rejection  of  all  lots.  Within  this  randomization  process  the  two 
point  OC  definition  problem  Is  cast  as  a  linear  program.  Employing  the 
traditional  solution  and  the  LP  solution,  a  procedure  is  devised  for  an 
exact  solution  to  the  problem.  ^ 
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1.  INTRODUCTION 


A  well  known  problem  In  lot  by  lot  attrloure  sampling  is  the  deter¬ 
mination  of  a  single  sampling  plan  whose  OC  curve  passes  through  two 
predetermined  points.  Usually  the  points  selected  are  the  process  averages 
and  the  corresponding  a  and  8  values  representing  levels  of  the  producer's 
and  the  consumer's  risks.  Ordinarily  this  problem  cannot  be  solved  exactly 
as  the  number  couple,  (n,c),  defining  the  sampling  plan  must  consist  of 
Integers.  An  easily  acquired  approximate  solution  to  the  problem  and 
tables  for  its  implementation  are  given  In  Grubbs  C33  and  Duncan  121. 

A  reformulation  of  the  problem,  which  Is  soluble,  requires  the  deter¬ 
mination  of  (n.c)  such  that  Pa(pj)  >  1-a  end  Pa(p2)  <0  for  the  smallest 
possible  c.  The  function  Pa(p)  defines  the  OC  curve  of  the  plan  and  p^ 
and  p2  ere  values  of  the  process  average.  It  is  assumed  that,  p^<  p? 
and  1-a  >  e.  Hald  14,53  offers  a  lucid  survey  of  both  exact  and  approxi¬ 
mate  solutions  to  this  problem.  Stephens  E73  derives  a  closed  form  solu¬ 
tion  to  the  exact  problem  by  employing  a  normal  curve  approximation  to 
the  binomial  distribution  given  by  Bbrges  [13.  The  solution  of  Stephens 
Is  extended  by  Jaech  C63. 

The  solutions  referred  to  represent  a  compromise  since  Ideally  we 
seek  a  sampling  plan  whose  OC  curve  satisfies  the  two  points  exactly. 

In  this  paper  we  associate  pre-determlned  probabilities  with  each 
outcome  of  the  sampling  process.  The  sampled  lot  Is  then  accepted  or  re¬ 
jected  accordingly  as  a  random  event  occurs  or  falls  to  occur  based  on  these 
probabilities.  Ne  refer  to  this  procedure  as  the  randomized  acceptance  plan 
Within  this  randomization  procedure  we  can  determine  the  probabilities 
which  cause  the  OC  curve  of  the  plan  to  pass  through  any  two  points  exactly. 


2.  RANDOMIZED  ACCEPTANCE  SAMPLING  PLAN 

Ue  propose  to  develop  a  randomized  acceptance  sampling  plan  In  tdtlch 
a  given  lot  Is  accepted  or  rejected  accordingly  as  a  random  event  occurs 
or  falls  to  occur,  respectively.  To  this  end,  consider  the  set  of  manbers 
{  *q.  ...»  *n }  •  where  0  ^  ^  - 1  for  1  *  0*  1*  ....  «-1. 

Subsequent  to  extracting  a  random  sample  of  size  n  from  the  lot.  the  ran- 
domlzed  plan  Is  executed  by  observing  x.  the  number  of  defectives  In  the 
sample,  and  accepting  the  lot  with  probability! or  rejecting  the  lot 
with  probability  t  .  Just  as  the  traditional  OC  curve  defines  the  average 
fraction  of  accepted  lots  for  a  given  level  of  the  process  average,  the 
mathematical  expectation  of  (1-t  )  defines  the  average  fraction  of  lots 
which  are  accepted  under  the  randomized  plan  for  a  given  level  of  the  pro¬ 
cess  average.  Therefore,  while  Pt  defines  the  operating  characteristics 
of  the  traditional  plan,  the  value  E(1-y  defines  the  operating  charac¬ 
teristics  of  the  randomized  sampling  plan. 

The  traditional  plan  Is  defined  by  (n,c)  and  the  randomized  sampling 
plan  Is  defined  byCVg,  tj,  *n  }  . 

In  the  notation  of  the  randomized  plan  the  traditional  plan  becomes 

{0,  for  x  •  0,  1,  ...»  c 
1,  for  x  »  c*0,  ....  n. 

Obviously,  under  the  randomized  plan  as  defined  above, 

so-y  •  o-y  *0  ♦ ...  ♦  o-y 

4  b,  4  ...4bc 

•v 

where  bx.  x  *  0,1,..., n,  Is  the  probability  mss  function  of  x,  the  number 
of  defectives  In  the  sample  of  size  n.  Usually  bx  will  be  defined  by  the 
binominal  or  Poisson  model. 


,iv| *  <  1 


3.  THE  TWO  POINT  OC  PROBLEM  AS  A  LINEAR  PROGRAM 

• 

Within  the  randomized  sampling  plan  concept  we  define  the  two  point 
OC  curve  problem  as  follows.  Having  selected  the  points  (p^,  l-o)  and 
(p2»  *)•  Pj«  P2.  1*D>  ••  **nd  •  *et  of  numbers  Pq,  ...»  *n)  such  that* 

‘W’-  (1) 

E2(1.tx)  •$ 

subject  to  0  <  Tj  <t1+1<  1,  1  •  0,1,...,  n-1.  The  symbols  Ej  and  E2 
denote  mathematical  expectations  calculated  at  p^  and  p2.  An  approximate 
statement  of  (1)  Is  obtained  by  casting  the  problem  as  a  linear  program. 


thus. 


subject  to: 


min  Sj  ♦  Vj  ♦  S2  ♦  V2 


Voi  Vn  ♦  •  •  •  ♦Vni  4  h  *  vi  *  •  \  (3) 

r0b02  Vl2  4  •  •  •  4V«2  *  S2  4  V2  ■  1  •  *  J 


V 


Tl* 


«  0 


1  *  0,  1,  ••••  n, 
where 


b  * 


n-1  *  \ 


Ti  *sr  V  S2*  V2 


« 


i  o. 


J  i  x*0»l *..»,••• •  • 


and  2j  ■  npj,  J*l,2. 

We  have  chosen  to  model  b^  as  Poisson  although  the  binomial,  the 
hyper geometric  or  any  relevant  probability  mass  function  on  x  would  be 
appropriate. 


The  variables  Sj  and  Sg  are  non-negative  error  variables  measuring 
the  extent  to  which  the  OC  curve  of  the  randomized  acceptance  plan  exceeds 
1-aat  pj  and  lies  below  a  at  respectively.  Analogously,  the  variables 
Vj  and  V£  ere  error  variables  which  define  the  distances  the  OC  curve 
of  the  randomized  plan  lies  below  1-a  at  pj  and  above  e  at  Pg.  It  Is  appa¬ 
rent  that  Sj  and  cannot  be  positive  together.  A  similar  statement  Is 
Implied  for  Sg  and  Vg. 

The  constraints  defined  In  (3)  accoranodate  OC  curve  departures  from 
the  equality  defined  In  (1)  where  the  constraints  defined  In  (4)  represent  a 
general  characterization  of  the  usual  properties  of  the  operating  characteris¬ 
tic  curve.  In  the  language  of  linear  programing,  the  constraints  of  (5) 
are  the  usual  IP  non-negativity  restrictions.  Setting  yq  •  0  and  vn  ■  1 
assure  the  OC  curves  of  the  randomized  plan  to  pass  through  (0,1)  and  (1,0) 
under  the  binomial  and  hypergeometric  model  assumptions.  All  these  constraints 
are  compatible  with  the  properties  of  the  traditional  OC  curve  and  are  In¬ 
tuitively  acceptable. 

Within  the  structure  of  this  problem,  our  goal  Is  more  extensive  than 
the  apparent  minimization  of  the  objective  function  as  stated  In  (2).  In 
addition  to  seeking  the  set{TQ,Yj,  •••»*„  I  for  which  $j  ♦  ♦  Sg  ♦  Vg  is 

least  we  also  desire  the  smallest  value  of  n  for  which  the  sum  of  the  errors 
is  equal  to  zero.  Surely  all  values  of  n  (n>0)  produce  basic  feasible  so¬ 
lutions  to  the  linear  programming  problem  as  stated,  so  we  turn  to  the  problem 
of  finding  {^q,  ...»  for  minimum  n  such  that  Sj  4  Vj  ♦  Sg  ♦  Vg  Is 

equal  to  zero. 

4.  DETERMINATION  OF  INITIAL  n 

From  the  Information  In  Grubbs' table  we  easily  obtain  four  traditional 
sampling  plans,  all  approximate  solutions  to  the  two-point  OC  curve  problem. 


Of  these  four  sampling  plans,  the  OC  curve  of  one  passes  through  (p^ ,  1 .  a) 
and  lies  above  a  at  p2.  Another  plan  has  an  OC  curve  which  passes  through 
(P2»  a)  and  lies  below  1-oat  p^.  We  designate  these  OC  curves  as  less 
discriminating  than  the  OC  curve  of  the  solution  plan  we  seek.  On  the 
other  hand,  of  the  two  remaining  approximate  solution  sampling  plans,  one 
has  an  OC  curve  which  passes  through  (pj,  1-w)  *nd  lies  below  8  at  p2, 
whereas  the  other  passes  through  (p2,  8)  and  lies  above  !-•  at  pr  We 
designate  the  OC  curves  of  these  two  sampling  plans  as  more  discriminating 
than  the  OC  curve  of  the  solution  sampling  plan  we  seek.  The  two  less  dis¬ 
criminating  OC  curves  are  characterized  by  sampling  plans  which  have  smaller 
sample  sizes  than  those  of  the  more  discriminating  OC  curves.  We  feel  In¬ 
tuitively  that  the  randomized  plan  we  seek  Is  defined  by  a  swnple  size  In¬ 
termediate  to  the  sample  sizes  of  the  two  groups  of  OC  curves  just  discussed. 
Therefore,  as  a  starting  point  or  an  Initial  sample  size,  we  select  the 
larger  of  the  two  sample  sizes  of  the  less  discriminating  OC  curves.  From 
this  Initial  n  we  will  Increment  upward  by  trial  and  error  until  we  reach 

the  least  n  for  which  Sj  ♦  ♦  $2  ♦  V2  Is  equal  to  zero. 

$ 

5.  THE  OPTIMUM  SOLUTION:  AN  EXAMPLE 

In  this  section  we  shall  Illustrate  the  solution  procedure  by  an 
example.  The  problem  we  have  selected  Is  presented  by  Duncan  (1074) 
and  Is  defined  as  follows.  We  seek  a  sampling  plan  which  will  reject  sub¬ 
mitted  material  which  Is  IS  defective  5%  of  the  time,  while  accepting  sub¬ 
mitted  material  tftlch  Is  K  defective  only  10S  of  the  time.  In  other  words, 
we  seek  a  sampling  plan  whose  X  curve  passes  through  the  two  points,  (0.01, 
0.95)  and  (0.06,  0.10).  Employing  Crubbs'  table  C33,  Duncan  C21 


defines  four  traditional  sampling  plans  which  approximate  the  desired  . 
plan.,  the  characteristics  of  these  four  plans  are  listed  In  Table  1, 
together  with  their  AOQL  values. 


Table  1.  Approximate  Solution  Sampling  Plans  Passing 
Through  (0.01 •  0.95)  and  (0.06,  0.10) 


Plan  no. 

n 

c 

n 

4 

■EEW 

1 

36 

1 

0.05 

0.23 

0.022 

2 

49 

1 

0.08 

0.10 

0.016 

3 

67 

2 

0.03 

0.10 

0.019 

4 

82 

2 

0.05 

0.04 

0.014 

*  Based  on  lot  size  of  1000 


It  Is  apparent  from  our  previous  discussion  of  section  4  that  the 
first  two  plans  of  Table  1  are  less  discriminating  than  the  plan  we  seek, 
while  plans  3  and  4  are  more  discriminating  than  the  desired  plan.  This 
means  we  Intuitively  seek  an  optimum  solution  with  a  value  of  n  Intermediate 
to  49  and  67. 

He  Initiate  a  solution  to  the  IP  problem  defined  by  (2),  (3),  (4)  and 
(5)  with «  ■  0.05,  B  *  0.10,  pj  ■  0.01,  p2  ■  0.08' and  an  Initial  value  of 

n  of  49.  LP  solutions  beyond  n  •  49  are  executed  by  Increasing  n  by 
one  over  the  previous  solution.  These  solutions  are  continued  until 
the  objective  function  equals  zero.  Table  2  contains  solution  data 
on  sample  size,  number  of  simplex  Iterations  for  each  solution  and 
the  values  of  each  error  variable  at  each  LP  solution. 
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Table  2.  LP  Results  for  the  Randomized  Plan 
(a  •  0.05*  6  *  0.10,  Pj  ■  0.01,  p2  •  0.08) 


Sample 

Size  No.  of  Simplex 

Solution  No.  n  Iterations  Sj  V1  S2  V2 


1 

49 

53 

— ywriVtl 

0.0361 

0.0000 

0.0000 

2 

50 

54 

0.0359 

3 

51 

55 

0.0355 

4 

52 

56 

0.0348 

5 

53 

57 

0.0338 

6 

54 

58 

0.0326 

7 

55 

59 

0.0310 

8 

56 

60. 

0.0291 

9 

57 

61 

0.0268 

10 

58 

62 

0.0242 

11 

59 

63 

0.0211 

0.0000 

12 

60 

64 

0.0 

EJ 

0.0167 

13 

61 

65 

0.0120 

14 

62 

66 

0.0073 

15 

63 

68 

0.0028 

16 

64 

68 

0.0000 

0.0000 

0.0030 

The  values  of  the  y’s  and  n  which  define  the  optimum  randomized  plan 

are: 


yQ  *0.000 
y,  *0.007 
y2  •  0.189 
y3  •  1.000 
• 

Y^  •  1.000 

n  •  64 


(6) 


l 


F 


V  v,V,*«V 


8- 


i 

r 


•  f 

t 

i 


It  Is  Apparent  that  sixteen  LP  solutions  were  required  to  arrive 
at  the  value  of  n  which  reduced  each  of  the  error  variables  to  zero. 

Tnls  number  was  arrived  at  by  Incrementing  upward  from  n  *  48.  It 
would  have  been  more  economical  In  practice  to  have  begun  with  n  ■  67 
and  Incremented  downward  to  the  optimum  n  at  64,  requiring  only  four  LP 
solutions  to  arrive  at  the  optimum.  However,  the  time  required  to  execute 
a  single  LP  solution  Is  relatively  small.  The  LP  computation  will  be 
discussed  In  section  7. 

The  randomized  plan  just  derived  Is  executed  by  observing  x,  the  num¬ 
ber  of  defectives  In  the  sample  of  size  64,  and  accepting  the  lot  If  x  *  0. 
If  x  *1,  a  biased  coin  which  falls  heads  with  probability  0.007  Is  tossed. 
If  the  coin  falls  heads  the  lot  Is  rejected,  otherwise  It  Is  accepted.  If 
x  *  2,  a  biased  coin  which  falls  heads  with  probability  0.189  Is  tossed. 

If  the  coin  falls  heads,  the  lot  Is  rejected,  otherwise  It  Is  accepted. 

For  all  values  of  x  equal  to  or  greater  than  three,  the  lot  Is  rejected. 
The  OC  curve  for  the  randomized  plar  Is  presented  In  Figure  1.  He  are 
assured  that  the  OC  curve  of  this  plan  passes  through  the  points  (0.0 1, 
0.95)  and  (0.08,  0.10),  exactly. 

The  similarity  between  the  randomized  plan  and  the  traditional  plan 
Is  apparent. 


6.  THE  AOQL  FOP  RANDOMIZED  SAMPLING  PLANS  UNDER  RECTIFICATION 

Under  total  rectification  of  rejected  lots,  the  AOQ  of  the  randomized 
sampling  plan  Is  defined  as. 


AOQ 


.  Jt 


p  [O-Vb0  ♦  O-Vh  *  ■ 

♦  n-VS)  • 


(7) 


-9- 


where 

0  <  1,  i*o,  1,2,  ....  ii-l  «nd  n  Is  the  sample  size  from 
a  lot  size  of  N.  Employing  the  Poisson  probability  distribution  defined 
as. 


where  z  »  np,  gives 

A0Q  *  Sr  e~2  lo-V* 4  ^*yi^2  4  o-y2)  -fr  ♦••• 4 

_n*0 

4  (1'yn)-nr3.  (8) 

Setting 

•■ffl  ■  0  from  (8),  and  collecting  like  powers  of  gives, 

0  *  0*V  4  <l4Yo  “  2V2  4  {l42Yl  •  3V"2T“4  ••• 

r  •  *n  (1-y Jz1*41 

♦  Iu"Vi  *  - -h—  •  w 

which  is  a  necessary  condition  for  a  maximum  AOQ  or  AOQL. 

Substituting  the  values  from  (6),  which  define  the  optimum  randomized 
plan.  Into  19)  and  solving  for  z  gives,  z*  *  2.182.  Substituting  z*  for 
z  In  the  expression  for  the  AOQ  defined  In  (8)  gives  the  AOQL  of  the  rando¬ 
mized  plan  as  0.0184.  It  Is  apparent  that  the  AOQL  of  the  randomized 
plan  Is  compatible  with  the  AOQL  values  of  the  approximate  solution  plans 
which  are  presented  In  Table  1.  The  AOQ  curve  of  the  randomized  plan  Is 
shown  In  Figure  2. 

Under  randomized  plans  Involving  rectification,  rejected  lots  are 
screened  1002  just  as  they  are  under  traditional  plans. 


Figure  2 •  AOO  Curve  -  Randomized  Plan 
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7.  COMPUTING  THE  LP  SOLUTIONS 

The  generation  of  the  linear  programing  model  for  computing  purposes 
was  done  by  a  matrix  generator  subroutine  written  for  UNDO  (Linear  Inter¬ 
active  Discrete  Optimizer),  a  generally  available  linear  programnlng  com¬ 
puter  package.  The  subroutine  requests  the  user  to  specify  n,  pj,  a,  p2 
ends,  and  from  these  Input  data  automatically  generates  the  desired 
linear  programming  model.  With  the  use  of  this  system,  generating  the 
randomized  sampling  plans  for  various  values  of  n  Is  done  Interactively, 
thus  model  Input  and  optimum  solution  require  only  a  few  minutes.  This 
subroutine  Is  available  from  the  authors. 

8.  CONCLUDING  REMARKS 

We  have  Introduced  the  notion  of  a  randomized  sampling  plan,  a  general 
concept  for  a  single  sampling  Inspection  procedure  by  attributes.  Within 
the  randomized  sampling  concept  we  have  cast  the  two  point  OC  curve  problem 
as  a  linear  program,  resulting  in  an  exact  solution  to  the  problem.  In 
addition,  we  have  demonstrated  how  to  determine  the  AOQL  for  a  specifically 
defined  randomized  sampling  plan.  Thus,  we  have  demonstrated  the  adapta¬ 
bility  ©f  the  randomized  plan  to  sampling  Inspection  with  rectification. 

In  addition  to  solving  the  two  point  OC  curve  problems  exactly,  the 
authors  feel  that  the  generality  and  ease  of  application  of  the  randomized 
plan,  together  with  the  power  of  the  LP  solution  technique  provide  a  useful 
and  ready  means  for  Investigating  many  relevant  problems  In  sampling  plan 
definition. 
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The  determination  of  an  operating-characteristic  (OC)  curve  which  passes  through 
two  arbitrarily  selected  points  is  a  well  known  problem  in  statistical  quality 
control.  Ordinarily  this  problem  cannot  be  solved  exactly  as  the  solution  — »pHi«g 
plan  must  be  defined  by  a  nonnegative  integer  couple  (n,c). 
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Block  #20  (continued) 

In  this  paper  a  randomisation  process  is  defined  far  the  acceptance  and  rejection 
of  all  lots.  Within  this  randomisation  process  the  two  point  OC  definition  problem 
is  cast  as  a  linear  program.  Employing  the  traditional  solution  and  the  LP  solution, 
a  procedure  is  devised  for  an  esact  solution  to  the  problem. 


